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Розглянуто основні обмеження задачі Алексідзе. На її підставі подано економічний стпо- 
сіб розв'язання задачі ітераційного уточнення трансформації магнітного поля, вигодя- 
чи з адитивного зображення магнітного потенціалу. Доведено теорему єдиності роз- 
в'язку вказаної задачі. 


У процесі вивчення особливостей трансформації регіональних та глобальних аномалій по- 
ля сили тяжіння (7) у зовнішньому просторі у" виявлено, що їх з адекватною потребам 
геофізичної практики точністю доцільно здійснювати шляхом розв'язання нелінійної гра- 
ничної задачі Алексідзе! (| для рівняння Лапласа, в граничних умовах якої безпосередньо 
стоять значення сили тяжіння 4(2) - |єгад И/ (2.)|. Узагальнена постановка задачі Алексі- 
дзе така: знайти потенціал Й/(т), х Є у", який задовольняє всередині замкнутої області 
у" з у" | )ду рівнянню Лапласа ДИ (т) - 0, х Є у", а в точках межі ду Ляпунова облас- 
ті у" і в нескінченності -- умовам: 
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де 4(х) - задана функція; у  - обмежена область точок евклідового простору В, яка 


зайнята масами Землі; у" -- її необмежене доповнення без мас, які тяжіють; ду -- межа 


областей у" й у" (поверхня вимірювань). 

Її розв'язок, у силу обраної однорідної шаруватої моделі середовища, визначається по- 
тенціалом простого шару з невідомою густиною, поширеною на контактній поверхні типу 
Ляпунова. Ця густина визначається як розв'язок рівняння сили тяжіння |?|, аналітичний 
вираз якого є наслідком аналізу властивостей модуля градієнта потенціалу | тай Й/(2))|, 
а конкретний вигляд залежить від обраної моделі Землі. У плоскому випадку ця задача 
отримала алгоритмічне розв'язання за допомогою ітерацій |3|. Точність її розв'язку істотно 
залежить від міри зумовленості цієї задачі, параметрів числового методу (числа і розта- 
шування фундаментальних розв'язків), наближень напряму шуканого градієнта, напряму 
зовнішньої нормалі, уздовж якої обчислюють похідні у граничних умовах (тобто врахування 
векторної природи поля) |4|. Однак її розв'язання успішне на простих моделях середовища 
за певних обмежень на поведінку розв'язку |4Ї. 

Задача (1) орієнтована на аналітичне продовження поля поза джерела аномальних мас 
при заданій з певною точністю моделі середовища, відомих обмеженнях на середовище та 
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чевидно, вперше задачу для поверхні /Іяпунова оприлюднено в доповіді А.Б. Чорного "задача Алек- 


сидзе для уравнения Лапласа и использование ее решений в геофизике / / Теория и практика интерпретации 
потенциальньтх полей", Ленінакан, 1986. 
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похибки вимірювань. На параметри середовища накладають обмеження у вигляді куско- 
вої неперервності і диференційовності за Фреше -- зображення приростів розв'язків задачі 
лінійною комбінацією приростів параметрів моделі, тобто виділення головної лінійної час- 
тини шуканого розв'язку с (Є). 

Основне апріорне припущення щодо поведінки шуканого розв'язку оберненої задачі -- 
лінійна незалеожність оператора, розв'язку та його похідних -- забезпечується одним з ме- 
тодів |5|. Указаних обмежень достатньо для забезпечення числової збіжності до точного 
розв'язку. Виявлена неоднозначність розв'язків задачі через їх скінченну вертикальну роз- 
дільну здатність, неадекватний вибір моделі задачі та/або початкових наближень розв'язку 
с(Є) властива самій природі обернених задач з похідними в граничних умовах (1). Точність 
дискретизації задачі залежить і від способу обчислень похідних. До нюансів числового 
моделювання належать: осциляція інтегральних ядер операторів оберненої задачі; істотна 
залежність візуалізації цих ядер від обраного способу гріддінгу. Нівелювати осциляцію мож- 
на, задіявши конструкції типу інтегралу Шварца для смуги |б|, а нюанси графічної візуалі- 
зації потребують експериментального дослідження та вибору адекватного інструментарію. 

Незважаючи на складнощі числової реалізації задачі Алексідзе, вже можна окреслити 
сферу Її застосування: уточнення фігури Землі та аналітичних трансформацій гравімаг- 
нітних полів. Зупинимось детальніше на способі відновлення магнітного потенціалу в задачі 
Алексідзе. 

Якщо магнітний потенціал зобразити як суму У(г) - (2) -- Т(2) нормального маг- 
нітного потенціалу М/(х) для даної епохи та збурювального магнітного потенціалу Т (2), то 
розв'язок задачі відновлення магнітного потенціалу У (2/)|г-о за заданими на границі зна- 
ченнями модуля його градієнта 7(2) - |зтад У (2), г Є ду, у спрощеному варіанті зведемо 
до такого ітераційного процесу відновлення збурювального потенціалу: 
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Далі визначаємо напрямні косинуси со5(п; 1, Ср) - у (2)/ іх) для визначення 
наближення модуля градієнта магнітного потенціалу: 
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який потрібен для обчислення чергового наближення збурювального магнітного потенціалу: 
Тіц2(2) З 222) сов(п1, т.) - (2) со8(м, т), х Є ду. (4) 
Звідси неважко знайти наближення магнітного потенціалу: 


УЧНІ (р) я (2) Туо(щ). (5) 
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За нормальний потенціал доцільно брати аномальне поле однорідно намагніченої сфе- 
ри; у такому разі збурювальний потенціал Т(2) відбиватиме поле реальних об'єктів, які 
не враховані у цій моделі. Крім даних (2), х Є ду (по-суті, це значення ДТ, (2)), нам 
відомі не лише рівняння фізичної поверхні ду Землі, а й вектор одиничної нормалі т(ї) - 
щ со5(хю, т.), К - 1,2,3 до цієї поверхні майже у будь-якій її точці. А оскільки відо- 
мо і нормальне поле (/ бо. ), то у будь-якій точці г замкнутої області у" можна обчисли- 
ти нормаль м(2) -- со5(м, є) до еквіпотенціальної поверхні (2) - Сх, яка проходить 

3 


1д 
(2) та значення со5(иу, т) - У; со8(м, а) со8(2к, т), де 
у дак Гаї 





через точку х, со5(и, тк) - 


1/ (2) з |втад У (2). 

Теорема 1. Якщо повержня Землі ду є множиною Ляпунова, то магнітний потен- 
ціал У(х), х Є у", даної спохи відновлюється однозначно з напруженості (2), х Є ду, 
магнітного толя цієї епоти. 

Схематично доведемо збіжність послідовних наближень У; (2), і - 0, со, до магнітного 
потенціалу реальної Землі У(2), х Є у". Припустимо, що вже збудоване і - 1 наближення 
як магнітного потенціалу Ус(2), У (27), . , Уг-1(2), так і напрямних косинусів со5(то, тр) - 
а сов(и, ту), оо, С08(пі, др), К - 1,3, 2 - 0,1,2, нормалі п(х) до еквіпотенціальної поверхні 
У(у) - Ст, яка проходить через розглядувану точку т Є 1". 

На подальшому і-му кроці послідовно визначимо 


з 
соз(та ті й со (т. р) со8(22;, т) 
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та граничну умову 
Б(2) - Х(2) сов(пі, т) - (2) со5(м, т) 


для знаходження і-го наближення густини д; (2) потенціалу простого шару Т; (2) з лінійного 
інтегрального рівняння Фредгольма: 
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Послідовність |7.(2)) збігається, що є наслідком із теореми, доведеної А.В. Чорним |?| 
(про природу потенціалу наразі не йдеться). 

Теорема 2. Якщо величиною є? (2) квадрату відношення модуля градієнта збурюваль- 
ного потенціалу до модуля нормального потенціалу можна знехтувати порівняно з є(2), 
то послідовність розв'язків 1Т.(2)) граничних задач 

АТ) об Є, ші Фо(2)-5Й, нлкоо В Ецц(є), сєду, (8) 
збігається до збурювального потенціалу Т(х), г Є у , якщо || 7(2.)|| с «Є || (2) с. 

Але виконання останньої умови при відновленні магнітного потенціалу зовсім не обов'яз- 
кове, оскільки гранична, умова 271(2)/дт., - Бці(2), х Є ду, є простішою, ніж гранична 
умова (2) - 2(2) сов(т;, т) -- (27) со8(м, т). Для збіжності послідовності 17(2)) достал- 
ньо, щоб існувала нерівність 2(2) «/ (2:)/ со5(и, т), хз Є ду. 

Зауважимо, що вказана схема відновлення потенціалу 7Т(х) значно економічніша за 
обсягом очислень, ніж схема, що запропонована в статті |ЗІ. 

Таким чином, достатньо, крім того, вказати, що якість числового розв'язку -- середньо- 
квадратичне відхилення (варіація) параметрів моделі задачі та вхідних даних від початко- 
вого наближення на поточній ітерації -- контролюється за критерієм нев'язки за додаткової 
умови дотримання мінімальної роздільної здатності. У будь-якому разі, вона залежить від 
міри зумовленості конкретної числової задачі та вибору початкових наближень розв'язку. 
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